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ABSTRACT 


To determine the feasibility of coupling the output of a single-mode optical fiber into a single-mode rib waveguide 
in a temperature varying environment, a theoretical calculation of the coupling efficiency between the two was inves- 
tigated. Due to the complex geometry of the rib guide, there is no analytical solution to the wave equation for the 
glided modes, thus, approximation and/or numerical techniques must be utilized to determine the field patterns of 
the guide. In this study, three solution methods were used for both the fiber and rib guide fields; the effective-index 
method (EIM), Marcatili’s approximation, and a Fourier method. These methods were utilized independently to 
calculate the electric field profile of each component at two temperatures, 20°C7 and 300°C r , representing a nominal 
and high temperature. 

Using the electric field profile calculated from each method, the theoretical coupling efficiency between an elliptical- 
core optical fiber and a rib waveguide was calculated using the overlap integral and the results were compared It 
was determined that a high coupling efficiency can be achieved when the two components are aligned The coupling 
efficiency was more sensitive to alignment offsets in the y direction than the x, due to the elliptical modal field profile 
of both components. Changes in the coupling efficiency over temperature were found to be minimal. 
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1 INTRODUCTION 


To model the behavior of an optical device over temperature, changes in both the refractive index and physical 
dimensions of the layers must be determined. Temperature-induced changes in the materials due to thermal expansion 
and stress are taken into account to determine device properties at elevated temperatures. 

This research explored the coupling efficiency between an elliptical-core optical fiber and a rib waveguide. The 
optical fiber transmits the source light into the guide, which serves as the input to an integrated optic pressure 
sensor. Because this sensor is to be used in a temperature varying environment on an airplane, changes in the coupling 
efficiency over temperature are of interest. The sensor design and fabrication were performed at the University of 
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Figure 1: Cross sections of rib waveguide and elliptical-core optical fiber. 


Cincinnati. 2 Both the rib and fiber physical parameters are illustrated in Figure 1. 


2 FIELD SOLUTION METHODS 


There is no exact solution for the field profile of a rib guide, however, there is an analytical solution for the fiber 
but it cannot be solved in closed form. 3 It was determined that the three techniques used to determine the rib field 
were also applicable to the fiber. 4 ~ 5 The effective-index method (EIM) and Marcatili’s method were chosen for the 
approximation techniques, and the Fourier method was chosen for the numerical technique, and was expected to 
produce the most accurate field results. All three methods provide solutions to the scalar wave equation, 


a 2 g x (x,y) 

dx 2 


+ 


d 2 gx(x,t/) 

dy 2 


+ [n 2 (x,y)& 2 -P 2 ] E x (x,y) = 0, 


( 1 ) 


where E x (x,y) is the field profile, n(x,y) is the spatially dependent function for the index of refraction, /? is the 

27T 

propagation constant, and k 0 equals — where \ 0 is the free space wavelength, 830 nm. The scalar-wave approximation 

K 

is valid if the difference between the guide core and cladding index is small, 


A = 


n j — 7X2 
fix 


« 1 , 


( 2 ) 


where ni is the core and ri 2 is the cladding index. 6 For the rib waveguide A = 0.064 and for the fiber A = 0.023, 
satisfying the weakly guiding condition. For quasi-TE modes propagating in the z direction, the electric field is of the 
form 


E (x,y,z,t) = x E x (x,y) e <(w ‘ 


(3) 


2.1 Effective-Index Method (EIM) 

2*1.1 Rib Waveguide - EIM 

The EIM 7-11 reduces the two-dimensional scalar wave equation, equation (1), into two one-dimensional problems. 
This is achieved by dividing the rib guide into separate planar guides as described in a previous paper. 12 Math- 
ematically, the EIM employs separation of variables to solve the scalar-wave equation, where the index profile is 
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approximated as 



n 2 (x,y)=n 

2 (x) + n" 2 (y). 

( 4 ) 

In equation (4), the quantity 

n 2 (y) is denoted as 




to 

II 

H-tO 

for y < |f> | 

( 5 ) 


n 2 (y) = n\ 

for y < — b 

(6) 


n ' 2 {y) = n\ 

for y>b* 

( 7 ) 

and n 2 (x) as 

n 2 (x) = 0 

for |x| < a 

(8) 


'2/ , \P 2 y 

71 (*) = - jfc|- 

n\ for |x| > a", 

( 9 ) 


where (3 y is the propagation constant solution for the y dependent field. Substituting the results for n 2 (x) and n ' 2 (y) 
into equation (4) yields the separable refractive index profile, n 2 (x,y), for which the EIM solves exactly. 13 The electric 
field profile of the rib is the product of the field distribution of the two one-dimensional planar slabs, X(x) and Y(v) 
expressed as v ’ Vi "’ 

E x (x,y) = X(x)Y(y). (10) 


2.1.2 Fiber - EIM 


The fiber shown in Figure 1 is a single-mode polarization-maintaining elliptical-core fiber. This fiber was approx- 
imated as a rectangle with dimensions a and b' such that the rectangle and ellipse have the same area and aspect 
ratio, given by 14 


a = a_ 2“ and 

Using the procedure given in section 2.1.1, the fiber field was 



determined for this rectangle. 


(ii) 


2.2 Marcatili’s Method 


2.2.1 Rib Waveguide - Marcatili’s Method 


For modes far from cutoff, the majority of the field energy is confined to the central core region. If this assumption 
applies, a pseudo-rectangular waveguide can be constructed to approximate the rib structure. 15 " 17 The cross section 
of the waveguide is divided into nine regions as shown in Figure 2, where t is the etch depth, r is the height of the 
etched SiON layer, such that the unetched height is r + t. 

Analogous to the EIM, this method reduces the two-dimensional wave equation into two separate one-dimensional 
problems as described in a previous publication. 12 The pseudo-guide 4is assigned a refractive index profile which 
approximates that of the actual rib and is defined by*^ 

n 2 {x,y)=n ,2 (x)+n" 2 (y)-n 2 , (12) 

where 


n'(x) = 
n'{x) — 


n i 
n 2 
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for |x| < W/2, 
for |x| > W/2, 


(13) 

(14) 




Figure 2: Marcatili’s method: pseudo-guide used to approximate rib guide. 


and 


for — t < y < r, (15) 

for ( y < -t) or ( y > r). (16) 

2.2.2 Fiber - Marcatili’s Method 


n"(y) = nj 

n"{y) = n 2 


Marcatili s technique can be used to calculate the field profile of an elliptical-core fiber. 4, ^ The elliptical core is 
approximated by a rectangle, for which an analytical solution can be found. The same rectangle dimensions chosen 
for the EIM are used here. 4 Using this rectangle as the pseudo-guide Marcatili’s method, as described in section 2.2.1, 
is used to determine the fiber field profile. 


2.3 Fourier Method 

2.3.1 Rib Waveguide - Fourier Method 


For the numerical technique, a solution to the scalar wave equation is obtained by expressing E x (x,y ) as an infinite 
series in terms of a complete set of orthonormal functions. 18-19 That is, 


Ex(z>y) — ^ ^ AprVpvfay), 


/i=l V=\ 


where 


^^{x,y) = sin (cr M z) sin (p v y ) , 


is a simple and complete orthonormal set of sine functions, and 


/i7T I/7T 

— — and p v — — , 

•L'x Lty 


(17) 


(18) 


(19) 
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where fi and v are integers. 
(0 < y < L y ), such that 


The sine function, is orthogonal over the finite base domain, (0 < x < L x ) and 


L , Ly 

II 

0 0 


^ nv{ x ,y)y ti'v'{x,y)dx dy = 


( 20 ) 


For this method, the field is only defined over this finite base domain; thus, proper choice of the domain limits is 
critical for accurate results. The values for L x and l/ y are chosen such that the field on those boundaries is zero. For 
well-guided modes these boundaries are small and this method performs well, however, for waves at or near cutoff 
this method does not apply. 


To obtain a matrix equation, equation (17) is inserted into the scalar wave equation, equation (1). Then, the left 
hand side is multiplied by tf M „(x,y) and integrated over the finite base domain producing 


m X TXly 

M=1 



( 21 ) 


where 


,/ii/ 


L X U 


L, Ly 

II 
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\ x >y) — 
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^(z,y)^v(z,2/) dx dy. 


( 22 ) 


In order to obtain a closed form solution for the integral, the index of refraction profile in equation (22) is approximated 
by a set of rectangles, each having a constant refractive index. 19 


Using the procedure from Marcuse, 1 ® two FORTRAN programs 29 were developed to determine the propagation 
characteristics of the rib and fiber. The L x and L y chosen values, 15 ym and 10 //m, respectively, were large compared 
to the fiber and rib guide geometry, but necessary to maintain the zero field condition on the domain boundary while 
the rib guide was shifted with respect to the fiber. This was done to determine the coupling efficiency at various 
transverse offsets. Increasing the number of terms, m x and m y , improves the accuracy until a convergence point 
is reached. Beyond that point, further increases in the number of expansion terms have negligible effect. It was 
determined that a series of 30 sine functions in both directions was sufficient to produce accurate results. 


2.3.2 Fiber - Fourier Method 


The Fourier method was applied to the elliptical-core fiber using the development in section 2.3.1 to calculate the 
field profile. The fiber core was approximated as a single rectangle with the integration taken over the elliptical-core 
boundary. The same L x and L y domain boundaries and expansion terms used for the rib were used for the fiber. 


3 Coupling Efficiency Calculation 

_ ^^ere are several causes of interconnect loss between a fiber and a waveguide. 21-23 These include mode field 
mismatch, transverse offset ( x f y offset), longitudinal separation ( z offset), angular offset, and refractive-index mis- 
match. The loss mechanisms addressed in this study were the transverse offset and mode field mismatch. Mode field 
mismatch contributes significantly to interconnect loss, 23-24 hence modeling is a necessary tool to design a waveguide 
whose field profile closely matches that of the input fiber. As misalignment in the x and y directions is inevitable in 
practice, the theoretical coupling efficiency results for x and y offsets are of interest. 

Once the electric fields of the rib guide and fiber were calculated, the coupling efficiency between the two compo- 


5 



ncnts, 77, is calculated using the overlap integral, 25 


*7 = 


J J Ef(x,y)E' g (x,y) dx dy 


J J E f (x,y)E’ f (x,y) dx dy J J E a (x,y)E*(x,y) dx dy 


(23) 


where the integrals are taken over all space and Ef(x t y) and E g {x,y) are the electric field profiles of the fiber and rib 
waveguide, respectively. For the EIM and Marcatili’s method, this integration is performed directly. For the Fourier 
method, the coupling efficiency is determined as the integral of the product of the fiber and rib guide fields squared, 
denoted 

(E E (24) 

Ap V and are the coefficients of the fiber and rib guide, given in equation (17), which are normalized as 

EE (25) 


The major axis of the elliptical-core fiber was oriented parallel to the x axis of the guide and the position of zero 
x and y offset was chosen to correspond to alignment of the maximum field positions of the fiber and rib. 26 After 
solving for the aligned case, various x offsets were introduced and a new coupling efficiency calculated for each offset 
at both temperatures. This procedure was then repeated for various y offsets while the x offset was fixed at zero. 


4 RESULTS 


4.1 Field Results 

Theoretical field results for the rib using the three methods are illustrated in Figures 3-5 with the rib cross-section 
superimposed. For all the contour plots in this paper, each contour indicates a region of equal field magnitude in 
increments of 0.1, with the maximum field magnitude being equal to 1.0. 

Field results using the EIM were very similar to those using the Fourier method, yet required considerably less 
computation time. The difference between the fields calculated using the EIM and Fourier method is shown in Figure 
6. Marcatili’s method assumes a tightly confined field, hence field results were tightly confined. The difference between 
the field calculated using Marcatili’s method and the Fourier method is illustrated in Figure 7. 

Modeling the elliptical fiber as a rectangle, field results using the EIM and Marcatili’s method are shown in Figures 
8-9, with the fiber cross-section superimposed. Field results from the Fourier numerical technique are shown in Figure 
10. Differences between the electric field profile calculated using the EIM and Marcatili’s method with respect to the 
Fourier method are shown in Figures 11 and 12. From these graphs, it appears that the field calculated using the 
EIM more closely matches that using the Fourier method than the field calculated using Marcatili’s method. 
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EWptk ol-Cor* Fiber (20 C) 



Figure 8: EIM: elliptical fiber field contours. 



Figure 11: Difference - Fourier and EIM fields. 


Marcotfii’s Method 

Elliptic at- Cor* Fiber (20 C) 



X (microns) 



Figure 9: Marcatili’s method: elliptical fiber field contours. Figure 12: Difference - Fourier and Marcatili fields. 



Figure 10: Fourier method: elliptical fiber field contours. 
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4.2 Physical Changes Due to Temperature - Results 


Compared to the 20® C data, the rib and fiber parameters changed slightly when calculated at 300°C, taking into 
account thermal expansion and stress with the results given in Table 1. It was apparent that the substrate stretched 
the film causing the rib width to increase while its height decreased. The indices of the rib and fiber all increased 
with increasing temperature Temperature-; induced change in the rib indices was two orders of magnitude greater 
than that due to the photoelastic effect. 


Field contours for both the rib and fiber at 300° were nearly identical to those calculated at 20°C' ) using each of 
the solution methods, but are not shown due to space constraints. 


PARAMETER 

20° C 

300° C 

W(300°C) 

2.700 pm 

2.702 pm 

^overlayer (300° C7) 

0.509 pm 

0.509 pm 

t-unetched (300° C7) 

0.516 pm 

0.516 pm 

^etched (300° {?) 

0.426 pm 

0.426 pm 

^overlayer (300° C) 

1.453 

1.456 

n rib (300°C) 

1.553 

1.556 

^isolation{ 300° C7) 

1.453 

1.456 


Table 1: Rib Guide Parameter Changes With Temperature 


PARAMETER 

20°C 

300° C 

a 

1.300 pm 

1.302 pm 

b 

0.650 pm 

0.651 pm 

ft core 

1.484 

1.487 

Tl cladding 

rT450 

1.453 


Table 2: Fiber Parameter Change With Temperature 
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4.3 Coupling Efficiency Results 


The x direction offset ranged from —3 fim to +3 fim in 0.10 fim increments. Coupling efficiency results between 
the fiber and rib guide, calculated using the three different methods, are shown in Figure 13 as a function of this 
x offset. It was evident from this graph that the EIM and the Fourier method predict similar results. Marcatili’s 
method predicts the fields of both components to be more tightly confined in the x direction, producing a coupling 
efficiency more sensitive to an x offset. 



- 3.0 - 2.0 - 1.0 0.0 1.0 2.0 3.0 

x OFFSET (microns) 

Room T (20 C) 


Figure 13: Calculated coupling efficiency vs. x offset. 


COUPLING EFFICIENCY vs. y OFFSET 



Figure 14: Calculated coupling efficiency vs. y offset. 

Similarly for the y direction, keeping the x alignment, the fiber position was varied from -3 fxm to -f 3 fim from 
the position of alignment and the theoretical coupling efficiency versus y offset is shown in Figure 14. The coupling 
efficiency versus y offset results were nearly identical for the three methods. From the graphs in Figures 13 and 14, 
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METHOD 


MAXIMUM CE (300°C) 

EIM 

85.03% 

85.18% 

Marcatili 

88.16% 

~ 88.15% 

Fourier 

81.99% 

81.98% 


Table 3: Maximum Coupling Efficiency 


it is apparent that alignment is more critical in the y direction, since the rib field calculated is more tightly confined 
in the y direction than the x. 

The coupling efficiency results versus x and y offsets at 300°C were nearly identical to the room temperature 
results. The maximum coupling efficiency from the three techniques at both temperatures is given in tabular form 
in Table 3. From this theoretical study, it was determined that the predicted changes in the coupling efficiency with 
temperature are negligible and would be difficult to detect in practice. 


5 Summary 


Results indicate that the EIM was more accurate than Marcatili’s approximation for both the fiber and rib 
■waveguide. While results from the Fourier method were more accurate, the amount of time required to determine 
the coupling efficiency for x and y offsets was several orders of magnitude greater than that for the approximation 
methods. 

These results indicate that the device is likely to perform properly in an environment with a temperature fluctuating 
between the two temperatures investigated. The calculated maximum coupling efficiency was higher than expected. 
Changes in the coupling efficiency with temperature were negligible because temperature induced variations in both 
the rib and fiber indices and dimensions were much less than expected. 

The coupling efficiency can be increased by varying the rib waveguide geometry. By reducing the rib width to 
2.4 fim and the etched height outside the rib to 0.35 /rm, the coupling efficiency increased to 88%, using the Fourier 
method. This was 6% greater than that calculated using the present design. Hence, modeling is beneficial in the 
design process to improve the theoretical coupling efficiency. 
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